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Abstract. A piienomenological model for the KN — > KS reaction is suggested. The model includes s 
and u channel exchanges by A, S, X'(1385), and yl(1520) and s channel exchanges by above-threshold 
hyperonic resonances. Explicit expression for the propagator for a particle with spin 7/2 is presented. 
High-mass and high-spin resonances play a significant role in the process. We deal with the whole set of 
existing experimental data on the cross sections and polarizations in the energy range from the threshold 
to 2.8 GeV in the center-of-mass system and reach a good agreement with experiments. Applications of 
the model to other elementary reactions of S production and to S hypernuclear spectroscopy are briefly 
discussed. 
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1 Introduction 

Whereas S hyperons were discovered soon after A and S 
hyperons, our knowledge in the S ^ —2 sector remains 
much more limited. S production reactions are typically 
more complicated than those for yl or Z" production, and 
their cross sections are small. Production of S' in a bi- 
nary process is possible with kaonic beams {KN — KS), 
which are relatively less intense, while other projectiles, 
such as pions, photons, electrons, protons, inevitably lead 
to three- or even four-body final states. That is why the 
S production dynamics is still poorly studied. 

Modern experimental facilities providing intense beams 
open new possibilities for systematic studies of S produc- 
tion. Photoproduction 7p — > K^ K^ has been investi- 
gated at JLab by CLAS collaboration [1,2]. Higher statis- 
tics will become available after completing of data anal- 
ysis of the new run gl2 [3,4]. The intense kaonic beam 
is expected at J-PARC where S production from nuclear 
targets is announced as the highest priority task [5,6]. 
The elementary reaction KN KS with formation of 
ground as well as resonant S states can be also studied 
at J-PARC [7]. PANDA collaboration proposed a study 
of the pp — SS reaction in view of the FAIR project [8, 
9]. Therefore, one may anticipate that much more empir- 
ical information on S ground state production as well as 
on the spectrum of S resonances will become available 
soon. Note that so far the spin and parity are known for 
the ground (octet), first excited (decuplet), and only one 
more resonant S state. Total number of observed S states 
is 11 (and not all of them are established with confidence) 
[10], which is much less than for = and S* = — 1, and 
than predicted by quark models. 
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The KN — > KS reaction was studied experimentally 
in the sixties and seventies [11-25] when the bubble cham- 
ber technique was popular. Later, the reaction left the 
focus of experimentalists, and only one measurement of 
the forward cross section at one energy was performed at 
KEK [26]. Few related theoretical papers {e.g., [27-30]) 
were published also in the past decades and dealt with 
small portions of the existing data. In these papers, some 
interesting features were pointed out, but no systematic 
theoretical picture has been ever presented. 

Later, closely related problems of SN interaction and 
S hypernuclear dynamics attracted some attention. Data 
on the KN — >■ KS reaction were discussed [31] in view of 
formation of S hypernuclei by the {K, K) reaction from 
nuclei. The {K^,K~^) reaction on nuclei was studied the- 
oretically {e.g., [31-34]). The forward cross section of the 
elementary K~p — > K~^S~ process needed in the calcula- 
tions was taken from empirical data. 

The SN interaction was analyzed theoretically within 
unitary meson-exchange models of baryon-baryon inter- 
actions extended to the full baryonic octet [35-37]. The 
SN potentials derived from those models can be discrim- 
inated, in principle, by the comparison with S hypernu- 
clear data. However, such data [31,38-40] are rather scarce 
and ambiguous now, therefore, the theoretical parame- 
ters (particularly, S coupling constants) relevant to the 
S — —2 sector, are not yet approved empirically. Here, 
the first attempt to extract the SN interaction cross sec- 
tion from analysis of the rescattering of S produced inside 
a nucleus [41] should be mentioned. 

Data on photoproduction reaction jp -> K^K^S^ [1] 
were analyzed theoretically by Nakayama et al. [42] in a 
diagram approach with coupling constants for S derived 
from flavor SU(3) symmetry. Probably, this is the first di- 
rect check of consistence of deduced S coupling constants 
with data. 
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Fig. 1. Diagrams for the s (a) and u (b) channel exchanges. 



In this paper, we present a systematical analysis of the 
existing data on the KN — )■ KE reaction at cm. energy 
y/s < 2.8 GeV in various charge channels. We develop a 
phenomenological model, which describes all the available 
data fairly well. Our analysis discloses some features of the 
E production dynamics and also shows restrictions caused 
by the limited amount of the existing data. Some prelimi- 
nary results have been published in conference proceedings 
[43]. The information extracted from the simplest binary 
reaction KN — >■ KE can be helpful for further studies of 
more complicated reactions like 7p K~^K^S~ . 

In sect. 2, we present our model. The experimental 
data are introduced in sect. 3. Our results are presented 
and discussed in sect. 4. We finish with conclusion and 
some outlook. 



2 Model 

We adopt the standard diagram approach. The simplest 
diagrams for our process are presented in fig. 1. We stress 
that the one-meson t channel exchange is absent since no 
(nonexotic) doiible-strangeness meson exists. So only two- 
meson exchange is possible in the t channel. Thus, we deal 
with the s channel [fig. 1(a)] and u channel [fig. 1(b)] ex- 
changes giving some comments on the t channel exchange 
in subsect. 4.1. 

Our treatment is restricted to the Born approximation, 
which means that neither channel coupling nor initial and 
final state interaction is taken into account. More com- 
prehensive approaches incorporating the above-mentioned 
effects and the unitary condition are known in the litera- 
ture. They are reviewed, for instance, in [44] for kaonic re- 
actions and in [45] for pion photoproduction. We use here 
the simpler approach considering that this is adequate for 
the first systematic study of reaction KN KS. 

We consider three charge channels of the reaction 



K-p K+S- 

K-p k°e:°, 



and 



K-n^K°S-, 



(1) 
(2) 

(3) 



for which experimental data are available. Due to isospin 
relations, the three corresponding amplitudes are not in- 
dependent: 



M{K-n^ K^S-) 
= M{K-p -)• K+S-) +M{K-p -s- K°S°) 



(4) 



Note that S exchange is allowed for all the reactions in 
the s as well as u channel while A exchanges is impossible 
in the s (u) channel for the third (second) reaction. 

Effective Lagrangian for B{l/2+)Y{l/2)K{0-) ver- 
tices is chosen with pseudoveetor couplings (Hereafter, no- 
tation and conventions of [46] are adopted): 



L = 



BYK 



B-,',[^^}Yd^K + h.c., 



(5) 



where the upper (lower) symbol in the curly brackets 
stands for positive (negative) parity of Y, respectively, 
other symbols have their usual meanings. Propagators for 
1/2 baryons have the standard form: 



S{q) 



A + M 



(6) 



For S(l/2+)r(3/2)iC(0-) couplings, the form adopted 
is also conventional: 



L 



^b\']y, 

TOtt I 75 J 



d'^K + h.c. 



(7) 



As for the 3/2 propagator, several different prescrip- 
tions are known [47-49]. Their fundamental advantages 
and drawbacks were discussed numerously {e.g., [50-52]). 
Treating our model as a phenomenological one, we do not 
enter into this discussion. For u channel exchanges as well 
as for s channel exchanges by particles below the reaction 
threshold, we employ the propagator, which is shown [49] 
to be the proper inverse of the Rarita-Schwinger equation: 



3iq^ - + iMr) 



M2 



M 



(8) 



Possible dependence of resonance widths on energy is ne- 
glected. 

For s channel diagrams with above-threshold interme- 
diate particles, we choose the propagators according to 
the simple criterion as follows. Let us consider an isolated 
s channel diagram [fig. 1(a)], where F is a spin-3/2 res- 
onance lying above the KE threshold and decaying into 
KE. We require that the cm. angular distribution of the 
products corresponds to the proper Legcndre polynomial 
[Pi(cos(6')) for 3/2+ and P2(cos(6l)) for 3/2" for the spin- 
independent part of the amplitude and the corresponding 
associated Legcndre polynomials for the spin-flip part] not 
only at the nominal resonance mass, but also at the slopes 
of the resonant peak. We checked that the only propaga- 
tor obeying this condition is one proposed by Adelseck et 
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al. [48]: 

^MQ) - s^q2_M^ + iMr) 
x(3,,.-7.>-^-^^^^^^^). (9) 

The same criterion is applied for higher spin resonances. 

A lot of Y resonances is known in the relevant en- 
ergy range. PDG compilation [10] gives eight resonances 
of four- and three-star status with 1.89 < M < 2.35 GeV. 
On the other hand, branching ratios of KS decay arc de- 
termined for none of them. Only for two of them, upper 
limits [3% for yl(2100) and 2% for ^(2030)] are deduced 
[10]. 

Probably, the branching ratios are small for all the 
resonances, while the main decay channels are ttA (for 
E resonances), ttE, and KN. This is physically natural 
since, first, the KE decay requires creation of an addi- 
tional ss pair and, second, the resonances lie not far from 
the KS channel threshold. On the other hand, cross sec- 
tions of the KN KB reaction are also not large, so 
even small branching ratios can contribute sizably to the 
reaction. Thus, the role of the above-threshold resonances 
should be studied. Most of these resonances has high spins. 
Therefore, we need technicalities for 5/2 and 7/2 spins of 
exchanged particles. 

For a 5/2 spin, we use the effective Lagrangian 

i= Z£22£Bpny a/^a^/^ + Zi.c. (10) 
L 1 J 

and the propagator from [53]: 

a ^ A + Vs 

10(g2 - M2 + iMT) 

~l~ PlXpl^l P<TIJ,'Pvi'' ~^ Pvp'y^l P<TV' Pp,p.' 

PaiJb' Pfj,!/' m) 

which obeys the criterion formulated above. Here P^i, = 

-9liv + <liJ.<lu/s. 

We construct the effective Lagrangian for a 7/2 spin 
similarly to those for lower spins: 

L = i^^B I ^ 1 Y^,pd''d''dPK + h.c. (12) 

An explicit expression for the 7/2 propagator were pro- 
posed recently [54]. We obtained and used another ex- 
pression, which is displayed and commented in Appendix. 

Apart from the ingredients discussed above, a form fac- 
tor suppressing high-momentum contributions is needed, 
which is especially important for high-spin exchanges. We 
employ the simple Gaussian form 

F(g) = exp(-qV^') (13) 

used in [55] within another approach. Here q is the cm. 
3-momentum of K ot K. This form factor is involved to 
each vertex in the diagrams of fig. 1. 



Form (13) differs from monopole, dipole, and other 
widely applied forms. We examined also other forms and 
found that form factor (13) is most successful for the data 
description. It reduces the cross sections at higher energies 
stronger than other parametrizations, which is compatible 
with the data. 

Evidently, F{q) is the function of s. Strictly speaking, 
this factor cannot be considered as a vertex form factor 
for the u channel exchange diagrams. We adopt the same 
factor for the u channel amplitudes too considering that 
form factors in such models actually are phenomenological 
ingredients rather than fundamental quantities. Similar 
practice has already been used in the literature within 
other frameworks {e.g., [55,56]). 



3 Input data 

Data on the integral and differential cross sections and 
polarizations were published [11-25] several decades ago. 
The integral cross sections are compiled in [57]. They have 
been displayed in various ways. Some of them have been 
tabulated, some of them have been shown on graphs, and 
others have been presented as parametrizations (Legendre 
polynomial expansions) . We accept only the cross sections 
presented directly (digitally or graphically) for our fit- 
ting procedure. As for Legendre polynomial parametriza- 
tions, they determine the cross sections (and polarization) 
with uncertainties, which magnitude is not always clear. 
In some cases they give even unphysical (negative) cross 
sections in limited angular ranges. Therefore, we do not 
deal with the Legendre parametrizations in the fitting pro- 
cedure. However, we check consistency of our results with 
the data presented in this form (sec the next section). 

Most of the data used has been presented with ex- 
perimental error bars. (Note that, reading the data and, 
especially, errors from graphs, we possibly introduce small 
additional uncertainties.) When errors have not been dis- 
played, we simply estimate the statistical errors as VN, 
where N is the number of events. 

To fit our model parameters, we use the data as follows: 

— integral cross sections of reaction (1), 56 points [11-13, 
15,17-25]; 

— differential cross section of reaction (1) at 11 values 
of cm. energy (1.95, 1.97, 2.07, 2.11, 2.14, 2.24, 2.28, 
2.33, 2.42, 2.48, and 2.79 GeV), 234 points [15,18,20, 

21,23,24]; 

— integral cross sections of reaction (2), 30 points [11-13, 
16,18,25]; 

— differential cross section of reaction (2) at 8 values of 
cm. energy (1.97, 2.02, 2.07, 2.11, 2.14, 2.15, 2.28, and 
2.47 GeV), 76 points [12,15,16,18]; 

totally 396 points. 

Also a very limited sample of data on reaction (3) is 
available [12,14,25]. Since they influence the fitting results 
rather weakly we do not include them into the fitting pro- 
cedure. Instead, we use these data for independent check- 
ing of our model. 
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Table 1. Results of the fits described in the text: values, products of coupling constants (14), and cutoff parameters {A\ for 
the low-mass particles and A2 for high-mass ones). 





x' 


Fa 


Fs 


-Fi:(1386) 


-Fa(1620) 


yli (MeV) 


-^£(2030) 


-^£(2250) 


yl2 (MeV) 


Without high-mass resonances 


1121 


0.4502 


0.1766 


0.0340 


-0.6462 


773 








With high-mass resonances 


985 


0.3303 


0.1185 


-0.0057 


-0.4076 


839 


0.0203 


-0.0838 


440 



Similarly, we do not include existing data on S po- 
larization [18,23] to the fit. The rather crude polarization 
data also have almost no effect on the fitting results. We 
discuss the polarization data and their description in our 
model in subsect. 4.2. 



4 Results and discussion 

4.1 Fit with low-mass exchange particles 

In the first stage, we restrict ourselves to four low-mass 
exchange hyperons in the s and u channels. In this stage, 
we have five free parameters: four products of the coupling 
constants 

Fy = /nyk/syk (14) 

and one universal cutoff parameter Ai. All these parame- 
ters are determined by the direct fitting to the data. 

Wc performed the fit with four exchange particles: 
yl(l/2+), S{l/2+), i;(1385,3/2+), and /1(1520, 3/2") in 
the s and u channels. All these particles lie well below the 
threshold (1812-1819 McV depending on the charge chan- 
nel). Comparison of the model and experimental data is 
presented in figs. 2 and 3 for reactions (1) and (2), respec- 
tively. Results of this fit are shown in the figures by dot- 
dashed lines. The parameters obtained and the value 
are shown in table 1. 

We tried to vary the number of low-lying exchange 
particles. First, even the simple version with only two ex- 
change hyperons {A and E) reproduces the data quali- 
tatively. However, the quantitative agreement is poorer 
(X^ = 1683), particularly, the backward peak in the dif- 
ferential cross sections is strongly and systematically un- 
derestimated. 

On the other hand, when we add the fifth exchange 
particle, yl(1405), to the four ones, we do not get a mean- 
ingful improvement. Hence, we conclude that the four- 
particle scheme is practically optimal as the starting point 
of our model. Of course, it should be noted that the model 
is constructed on purely phenomenological grounds. We do 
not think that exchanges by other particles are truly ab- 
sent. Instead, the four exchanges reflect effectively a more 
complex real physical picture. 

It is seen from figs. 2 and 3 that the fit reproduces 
well the main features of the cross sections. The energy 
dependence of the integral cross sections is described well 
within the whole energy range with a physically reason- 
able value of cutofi^ parameter Ai. The backward peak in 
the differential cross sections is fairly reproduced though 
some underestimation at certain energies should be also 
noticed. Shape evolution of the differential cross sections 
with energy is also reproduced. 



We checked a possible contribution of the t channel 
exchange mechanism adding the diagram with one-boson 

exchange by a double-strangeness meson. This meson, of 
course, need not be assigned to any non-qq exotics, but 
rather may be treated as a simulation of two-kaon ex- 
changes. Varying reasonably properties of this effective 
meson, we did not obtain essential improvement. So we 
suppose that the t channel exchange plays a minor role. 

Obtained value x^ = 1121 at 391 degrees of freedom is 
not low. In our opinion, an important contribution to this 
vahie originates from experimental pitfalls. Compare, for 
instance, the experimental cross sections at energies very 
close to each other, ^/s = 2.14 and 2.15 GeV, from differ- 
ent experiments [15] and [16], respectively, in fig. 3. It is 
seen that they differ from each other strongly. Physically, 
such behavior can be explained only if a very narrow reso- 
nance exists in this region, which is rather unlikely. Prob- 
ably, no realistic theory can reproduce these cross sections 
simultaneously. Thus, one can suggest that the data are 
uncertain. Some other similar though less prominent self- 
contradictions of the data can be; indicated. So we think 
that the x^ value can not be reduced radically. Of course, 
this restricts possibilities to deduce a fully reliable theo- 
retical picture. 

On the other hand, some drawbacks of our descrip- 
tion appearing systematically should be noted. The in- 
tegral cross sections in the bump region (approximately 
1.9-2.1 GeV) are underestimated. We can not describe 
the forward peak, which is seen in reaction (2) [but not 
in reaction (1)]. Unfortunately, this peak is concentrated 
in the single bin at each energy, so its structure remains 
unclear. Possibly, our model in the four-particle version 
can not describe also the data on S polarization (see the 
next subsection). We tried various modifications of the 
model restricted to sub-threshold exchange particles. Par- 
ticularly, we examined the form factor depending on u in 
the u channel diagrams, which can affect the angular dis- 
tribution. We employed the form from ref. [42]. However, 
in this and some other ways we could not achieve essential 
improvement. So, we suggest that feasibility of the model 
with sub-threshold exchange particles is limited. An ex- 
tension of the model is presented in the next subsection. 

4.2 Incorporation of above-threshold resonances 

As mentioned above, eight well established Y resonances 
in the relevant energy range are known (and also more 
than ten doubtful ones). Of course, it is impracticable to 
incorporate all of them in the calculation since the avail- 
able data do not give a possibility to determine eight addi- 
tional parameters (products of coupling constants). More- 
over, masses and, especially, widths of the resonances are 
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Fig. 2. Integral (top panel) and diflterential (other panels) cross sections of the reaction (1). The solid lines are the results from 
the full version of the model. The cross sections calculated only with the s channel exchanges by X'(2030) and 17(2250) are 
shown by dashed lines. At ^/s > 2.4 GeV, the resonance contribution becomes negligible and dashed lines cease to be visible. 
The dot-dashed lines are obtained from the fit without higher-mass resonances. The experimental differential cross sections are 
taken from [24] (1.95 GeV), [15] (1.97, 2.07, and 2.14 GeV), [18] (2.11, 2.28, 2.42, and 2.48 GeV), [23] (2.24 GeV), [21] (2.33 
GeV), and [20] (2.79 GeV). 
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Fig. 3. The same as in fig. 2 for reaction (2). The experimental differential cross sections are taken from [15] (1.97, 2.07, and 
2.14 GeV), [12] (2.02 GeV), [18] (2.11, 2.28, and 2.47 GeV), and [16] (2.15 GeV). 



also uncertain and, strictly speaking, may be free param- 
eters too. 

We examined various combinations of several reso- 
nances from the eight known ones adding to the four- 
particle scheme as above and obtained in several cases sub- 
stantial improvement of the data description. Ultimately, 
we choose 17(2030) and 17(2250) stressing that this choice 
is not unique. Spin and parity (J'^ = 7/2+) of r(2030) 
are well established. For 17(2250), most probable assign- 
ments are 5/2" and 9/2" [10]. We choose J"" = 5/2" for 
17(2250) checking that other choices do not change the re- 
sults drastically. We fix the "nominal" masses (2030 and 
2250 MeV) and widths 175 and 105 MeV, respectively. 



The latter quantities are the median values of the ranges 
presented for the widths in the PDG compilation [10]. 

The corresponding form factor is chosen in the same 
form (13). However, we need to employ another cutoff pa- 
rameter A2 for the above-threshold resonances. The reason 
is as follows. The fitting requires for the low-mass parti- 
cles the values of cutoff parameter Ai about 1 GeV. If we 
calculate an isolated high-mass and high-spin resonance 
with such values of Ai, we obtain a peak shifted to the 
right by several hundreds MeV from the nominal mass. 
We consider such behavior as inconsistent with empirical 
data because the resonances are revealed as peaks at near- 
nominal masses. More hard cutoffs {A2 about 0.5 GeV) 
put peaks to their proper places. 
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Fig. 4. The same as in fig. 2 for reaction (3). The experimental differential cross sections are taken from [12] (2.02 GeV) and 
[25] (2.62 GeV). 



Results of our final fit are shown in figs. 2 and 3 by 
solid lines. The best fit parameters are presented in table 
1. It is seen that the description of the data is improved 
substantially though decrease is not so large. Particu- 
larly, the drawbacks mentioned in the previous subsection 
are remedied. First, a strong forward peak appears for re- 
action (2) in agreement with the experiments. The bump 
in the integral cross sections is described adequately now. 

The data manifest prominent forward peaks at i/s = 
2.0 — 2.3 GeV in the differential cross sections of reaction 
(2), but not (1). Our analysis gives a simple explanation 
of this feature. It appears that reaction (2) in the bump 
region is almost purely resonant (see dashed lines in fig. 3). 
The higher-mass resonances give nearly equal contribution 
to the cross sections of both reactions. However, reaction 
(1) is also substantially contributed by A and yl(1520) u 
channel exchanges, which are impossible for reaction (2). 
That is why the differential cross section of reaction (2) is 
more symmetric than that of reaction (1) and the "bump" 
integral cross section of reaction (1) is about twice as large 
as that of reaction (2). 

The integral and differential cross sections of reactions 
(1) and (2) considered as most reliable are used for the 
fitting. Besides, other data mentioned in sect. 3 exist and 
they are not included to the fitting. It should be noted that 
the total statistical weight of those data is relatively small 
so they anyway cannot affect the fitting strongly, consider- 
ably increasing amount of the calculations in some cases. 
We prefer to use these data for independent checking of 
our model. 



In fig. 4, the integral and differential cross sections of 
reaction (3) are displayed. It is seen that the addition of 
the higher-mass resonances improves the description con- 
siderably though accurate quantitative agreement is not 
reached. It is notable that just S resonances are needed 
for this improvement while A resonances in the s channel 
do not contribute to reaction (3). 

Then, we show in fig. 5 the differential cross sections 
restored from the Legendre polynomial expansions com- 
pared to our calculations. Since the uncertainties of these 
data are not well defined, the agreement can be estimated 
only visually. It is seen that the theory is reasonably con- 
sistent with the experiment. Considerable discrepancies at 
the near-threshold energies are possibly associated with 
final-state interactions neglected in our approach. 

Another sample of data is associated with S hyperon 
polarization. All the data that can be directly taken from 
papers [18,23] are presented in figs. 6 and 7. 

Without the above-threshold resonances, the polar- 
ization in our model is practically zero. Since we work 
within the Born approximation, nontrivial imaginarities 
in the amplitudes leading to nonzero polarization appear 
only from the widths in the propagator denominators. The 
subthreshold resonances have relatively small widths and 
can not give a significant polarization themselves. A siz- 
able polarization appears from the interference between 
the higher-mass and lower-mass exchanges. Note that at 
y/s > 2.4 GcV (beyond the resonance region) the polar- 
ization becomes almost zero again. 



D.A. Sharov et al.: Phenomenological model for the KN — > reaction 




--e- -e. 0.0 0.5 --e, -ei 0.0 0.5 --e, -e. 0.0 0.5 1.0 

cos^ 

Fig. 5. Comparison of our calculations (solid lines) with the experimental differential cross sections of the reaction (1) restored 
from the Legendre expansions [15] (dashed) and [12] (dot-dashed). 
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Fig. 6. Polarization of S hyperons times the differential cross 
section from reaction (1). Convention for curves is the same as 
in fig. 2. The experimental data are taken from [23] (2.24 GeV) 
and [18] (the other energies). 



Fig. 7. The same as in fig. 6 for reaction (2) . The experimental 
data are taken from [18]. 



The data are crude and even the hypothesis of exact 
zero polarization cannot be rejected. However, it seems 
that the final version of the model gives a better corre- 
spondence with the data. We consider this as one more 
argument for the important role of the higher-mass reso- 
nances though it is possible that effects neglected in our 
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Fig. 8. Polarization of S hyperons times the differential cross section from reaction (1) restored from the Legendre expansions 
[15] (dashed) compared with our calculations (solid). 



model (like initial and final state interactions) can change 
the picture substantially. 

In fig. 8, the polarization data restored from the poly- 
nomial expansions [15] are presented together with our 
calculations. Uncertainties of the data are not well defined 
again, so only visual comparison is possible. It is difhcult 
to deduce robust conclusions from this comparison. Note 
that the data are often in contradiction not only with the 
calculations, but also with the data at a very close energy 
(compare, for instance, the data at ^/s = 1.99 and 2.01, 
or 2.15 and 2.17 GeV). The largest discrepancies are seen 
in the region of backward angles where the polarization 
varies very sharply. (Note that the polarization at = 0° 
and 180° is, of course, exact zero, while the curves are 
sometimes almost vertical in these regions so it may seem 
visually that the polarization has some nonzero value.) 

As mentioned above, we tested also combinations 
of the higher-mass resonances other than 17(2030) and 
17(2250). Here we explain our final choice. 

Underestimation of the integral cross sections in the 
bump region pointed out in subsect. 4.1 requires a reso- 
nance with the mass about 2.0-2.1 GeV. We observed that 
yl(2100), yl(2110), and i7(2030) give similar quality of de- 
scription of the integral cross sections of the reactions (1) 
and (2). Our choice of 17(2030) is inspired by the reaction 



(3) (see fig. 4), where only E resonances can contribute 
in the s channel. It should be noted, however, that there 
is only one experimental point in the bump region for the 
integral cross section of this reaction. 

The significant role of the i7(2030) resonance in the 
KN — ^ KS reaction was suggested by us in [43]. It is 
remarkable that very recently the similar conclusion was 
drawn [58] for S photoproduction 7p — > E!~ . 

The second resonance, 17(2250), provides a somewhat 
better description of the cross sections at the right slope of 
the bump. Also substantially nonzero polarization appears 
in a wider energy range. Since the spin of 17(2250) is not 
yet established, we checked the sensitivity of the results 
to the spin value. Changing the spin to 7/2 and 9/2, we 
found only small modifications of the results, though the 
fit became a bit worse. 

We tried also to incorporate larger number of the 
higher-mass resonances (up to seven). Of course, the fit 
becomes more ambiguous and sometimes multiple min- 
ima of appear. More important, we did not find great 
benefits in the data description. So we conclude that the 
suggested model is adequate to the existing amount of the 
data. 

We calculated the partial widths of £'(2030) and 
17(2250) decays into KS using our products of the vertex 
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constants and known branching ratios for the KN chan- 
nel from [10]. We obtained 1.2 and 0.3 MeV, respectively. 
Though the total widths are known with large ambigui- 
ties, we can roughly estimate the corresponding branching 
ratios as 0.7% and 0.3%, respectively, which is not in con- 
tradiction with small upper limits deduced in [10]. It is 
remarkable that so small branching ratios imply in the 
cross sections significantly. 

Many authors tried to extract vertex constants for 
couplings NKA and NKS (for a recent compilation, see 
[59]). Some predictions for the SKA and SKS couplings 
are made within general SU{3) schemes {e.g., [35]) or in 
more specialized approaches [60,61]. Some fragmentary 
attempts to derive the constants for resonances are also 
known {e.g., [55,42]). 

Even for the NKA and NKS couplings, the values 
deduced by various ways are several times different [59]. 
Often only the coupling constants are presented rather 
than the full vertex functions while the coupling constants 
practically have a physical meaning only together with 
the relative form factor. For this reason, the comparison 
of our constants with highly ambiguous and widely scat- 
tered corresponding quantities from the literature is not 
instructive. Note that we determine only products of the 
coupling constants, and our form factors are not normal- 
ized to unity on the mass shell as seen from (13). 

We examined the set of the coupling constants from 
[42]. Incorporating also form factors from [42] and treat- 
ing the cutoff parameter as a free one, we can not achieve 
quantitative agreement with the data. The salient dif- 
ference between the two models is considerably stronger 
couplings -Fi;(i385) and ^a(1520) iii [42] with respect to 
our ones. To compare properly the couplings, we fixed 
the constants from [42] and fit parameters A of our form 
factor (13) allowing the cutoff parameters to be differ- 
ent for different hyperons. We obtain rather hard cutoffs 
{A = 400 - 500 MeV) for 17(1385) and yl(1520) contrary 
to the ground state A and S. This means that the data 
on the reaction KN — )■ KS require to reduce these vertex 
functions. Note that our constants are purely phenomeno- 
logical while the constants in [42] are determined from 
SU (3) symmetry relations. 

Reaction {K~,K^) on nuclei becomes a hot topic 
now owing to the study of S hypernuclei starting at 
J- PARC [5,6]. Predictions of the cross sections of the 
'^7j{K~ , K^)^{'L—2) reaction {e.g., [32]) incorporate pro- 
ton effective numbers, which are determined by the nu- 
clear and hypernuclear structure, and properties of ele- 
mentary process (1). 

Generally, the amplitude of the elementary process has 
the form 



f{e)^g{e) + h{e)cT-n. 



(15) 



To our knowledge, all calculations for S hypernuclear pro- 
duction so far neglected spin-fiip component h{9) (with 
the single exception of ref. [62] where the h amplitude is 
chosen purely arbitrarily). In this case, the single input 
related to the elementary process is the differential cross 
section, which can be taken [32] directly from the data. 
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Fig. 9. Spin-independent and spin-flip amplitudes at kiab = 
1.8 GeV/c (-^s = 2.14 GeV) in the laboratory system at small 
angles. 



The spin-fiip component can lead to some important 
effects in hypernuclear production, for instance, to abnor- 
mal parity transitions as was considered earlier [63] for 
A hypernuclei. In view of hypernuclear applications, we 
present in fig. 9 amplitudes g and h in the laboratory 
system at incident momentum 1.8 GeV/c and forward an- 
gles. It is seen that spin-flip amplitude h grows quickly 
with the angle and becomes comparable with g even at 
Slab ~ 5°. Therefore, the spin-flip amplitude is generally 
non-negligible. Probably, g plays nevertheless a leading 
role, but a study of production of S hypernuclei with tak- 
ing into account spin-flip effects is worthwhile. 



5 Conclusion and outlook 

We have formulated the phenomenological model for the 
KN — KS process. While binary reactions for produc- 
tion of A and S hyperons are studied thoroughly, this 
process was completely out of theoretical scope for the 
last decades. Though the quality of the data is far from 
adequate, dealing with the whole set provides a possibility 
of the quantitative analysis of the reaction. 

Our goal was to elaborate a working tool for treating 
the reaction of S production rather than to consider it 
from the fundamental point of view. Therefore, the model 
should be regarded as effective one. We tried to minimize 
the number of exchange particles so the coupling constants 
possibly include effectively contributions of other hyper- 
ons left aside the model. Nevertheless, we hope that the 
model can be a sizable step for developing the theory of 
more complicated reactions of S production, particularly, 
photoproduction 7p — > K^K^S^ . The model suggested 
in [42] for the photoproduction fails to describe the data 
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on the KN — >■ KE reaction. This is not so surprizing be- 
cause the parameters of the model of ref. [42] as well as 
the parameters of our model are effective and probably in- 
corporate implicitly different neglected effects. Consistent 
simultaneous treatment of various reactions of E produc- 
tion is an ambitious future task. 

Earlier, ti channel exchange was considered as the dom- 
inant mechanism of the KN — >■ KE reaction. Confirming 
an important role of the u channel exchange, we show that 
the contribution of the on-shell resonances is also signif- 
icant. We cannot claim that just 17(2030) and Z'(2250) 
work (other choices are also possible), but establish that 
taking into account of the resonance mechanism is neces- 
sary. It is notable that relatively small (not higher than 
1%) branching ratios of y — > KE decays are manifested in 
the cross sections (which are not large themselves) clearly. 
We show that charge channel K~p — >■ K^E^ is "more res- 
onant" than K~p K^E~ due to lack of u channel A 
exchange in the former case. Maybe, in future this feature 
can be used for delicate studies of high-lying Y resonances 
in view of, for instance, J-PARC capabilities. 

Since the threshold of the reaction is relatively high, 
the majority of relevant resonances has large spins. Hence, 
we presented here technicalities for spin 7/2 exchanges, 
which can be helpful also for other reactions. 

As usual, polarization is the quantity rather sensitive 
to reaction mechanisms. Seemingly, our model catches the 
main trends of the polarization though more accurate data 
on polarization (as well as the cross sections) are needed 
for robust conclusions. 

The model is directly applicable to the theory of pro- 
duction of E hypernuclei. While empirical data on the ele- 
mentary process do not enable one to extract certainly the 
spin-flip amplitude, the model do. Predictions of hypernu- 
clear spectra measurable in the {K~,K^) reaction with 
account of the spin structure of the elementary amplitude 
and analysis of coming data from J-PARC is another task 
for the future. 



This work was supported in part by Russian Foundation for 
Basic Research, grant No. 08-02-00510. 



Appendix 



General recipes for propagators of particles with arbitrary 
spins were elaborated long ago [64-67]. While the propa- 
gator for spin 5/2 also has long been known [68], the 7/2 
case was treated explicitly [54] much later. Huang et al. 
[54] using the prescription of Behrends and Pronsdal [64, 
65] first derived the expression for the spin 4 projection 
operator -P,ii^2/i3^*4,i'ii^2i'3i'4 (4) and then got the 7/2 prop- 
agator as 



We repeated their derivation and obtained the follow- 
ing expression: 



(g 5Z "^''1 "'^ 1^2 PfiS 

{^lv} 

~84 S Pni'^iPn2'^2^i 



+ 



420 



•/is 1^3 



{ixu} 



Here 



and 



+ 



.(A.2) 



(A.3) 



(A.4) 



The sums in (A.2) are taken over all permutations of 

u indices in the first sum (6 terms) and over all permu- 
tations of /i as well as v indices in the other sums (36 
terms each). Due to simple identities Pafi — P^a and 
Pa^P-yS = P-ysPa^, the sccoud, third, and fourth sum con- 
tains 9, 18, and 9 nontrivial terms, respectively. 

The final expression from [54] can be simplified and 
reduced to the same form as (A.2). However, factor 
— 103/7560 appears in front of the second sum instead of 
our —1/84. We checked numerically that (A.l) and (A.2) 
(but not the formula from [54]) are equivalent. We also 
checked that propagator (A.2) satisficis the criterion from 
sect. 2 (gives correct angular distribution) while the prop- 
agator from [54] does not. So we believe that expression 
(A.2) is correct. 

Of course, one may use directly eq. (A.l). But in cum- 
bersome calculations, the explicit formula (A.2) is more 
convenient and time-saving. 

For the partial width of the channel ^(7/2=*=) 
K{0-)B{l/2+) we obtain: 



r 



YKB 



Eb±Me 



707rm6 My 



(A.5) 



where the sign in the numerator coincides with the par- 
ity of the resonance, q is the cm. 3-momentum of decay 
particles. Estimating the partial width in sect. 4, we mul- 
tiplied (A.5) by the relevant form factor and calculate q 
at the nominal mass of the resonance. 



q2_M-^+ iUr 9 



7 7 -fa;(LHjli2jli3,/3fll'2l^3 (4) 

(A.1) 
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